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Abstract 

The dynamical model on 3+1 dimensional space-time admitting soliton 
solutions is discussed. The proposal soliton is localized in the vicinity of a 
closed contour, which could be linked and/or knotted. The topological charge 
is Hopf invariant. Some applications in realistic physical systems are indicated. 
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1. Introduction 

The term "soliton" entered applied mathematics in 1965. It was coined by 
M. Kruskal and N. Zabusky for a special solution of nonlinear Korteweg-de Vries 
(KdV) equation, depicting solitary wave [jjj. Use of convention of particle physics 
language shows that the author envisioned the particle-like interpretation for the 
object which they called soliton. 

The attention of mathematical physicists to solitons was attracted after the 
inverse scattering method was devised by G. Gardner, J. Green, M. Kruskal and 
R. Miura for solving the KdV equation ||] and its extension to Nonlinear Schroedinger 
Equation was found by V. Zakharov and A. Shabat ||. In the 1970's this method 
and its generalizations got a lot of attention and involved quite a few active partic- 
ipants. Rather complete review can be found in (J]. In the end of that decade the 
quantum variant of the method was constructed and particle-like interpretations of 
solitons got natural confirmation in terms of quantum field theory, see review in || . 
The mathematical structure of the quantum method was deciphered in pure alge- 
braic way leading in the 1980's to notion of quantum groups with new applications 
in pure mathematics and mathematical physics. 

The value of solitons for the particle physics consists in the possibility of going 
beyond the paradigm of the perturbation theory. Indeed, soliton solutions corre- 
spond to full nonlinear equations and disappear in their linearized form. Charac- 
teristic for solitons is that they interact strongly if the excitations of the linearized 
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fields interact weakly. Another attractive feature is the appearance of elementary 
topological characteristics for solitons topological charges. 

This was understood already in the middle of 1970's by several groups as I 
underlined in my lectures, when I was touring USA in 1975 (see e. g. ||). However, 
all these tantalizing features of solitons had one very important drawback: the 
developed methods applied only in 1 + 1 dimensional space-time. 

Naturally the search for 3 + 1 dimensional generalizations became eminent. 
General considerations showed that many features of 1 + 1 dimensional systems, 
such as complete integrability and existence of exact many-particle solutions could 
not be generalized to 3 + 1 dimensions. However, the mere existence of "one- 
particle" soliton solutions was not excluded. One particular example was introduced 
by Skyrme in a pioneer paper long before the soliton rush. Another example was 
proposed by G. 't Hooft and A. Polyakov in 1975 0. In the following years their 
solutions got real applications in nuclear and high energy physics. 

In both examples the solitons are "point-like" , namely their deviation from the 
vacuum is concentrated around central point in space. Moreover they have spheri- 
cal symmetry, allowing the separation of variables in the corresponding equations, 
reducing them to ODE, which one can treat on a usual PC. 

In my lectures ||, already mentioned, I proposed one more possibility for 3 
+ 1 dimensional system, allowing solitons. The model, which superficially looks 
as a slight modification of Skyrme model, has quite distinct features. The center 
of the would-be soliton is not a point, but a closed contour, possibly linked or 
knotted. However my proposal remained unnoticed. The reason was evident: the 
maximal symmetry for such a soliton is axial, reducing 3-dimensional nonlinear PDE 
to 2-dimensional one. Existing computers were not able to treat such a problem. 
Thus my proposal was in slumber for 20 years until my colleague Antti Niemi 
became interested and agreed to sacrifice a year to learn computing and devising 
the programm. The preliminary results published in |^| attracted the attention 
of professionals in computational physics and now we have an ample evidence, 
confirming my proposal [Tc| ] , ]i"l| . 

The development which followed showed unexpected universality of my model. 
The variables, used in it, were shown to enter the list of degrees of freedom for several 
systems, having realistic physical applications |p^| , [fl3"| . 

In this talk I shall describe all these developments in detail. First I shall 
introduce the model, then briefly discuss its numerical treatment and finish with 
the description of the applications. 

2. The field configurations and Hopf invariant 

The space time is 4-dimensional Minkowski space M with linear coordinates 
a;' 1 , fj, = 0, 1, 2, 3, x° being time and x k , k — 1, 2, 3 space variables. The field n(x) 
is defined on M and has values on 2-dimensional sphere § 2 : 



n : M -► § 2 . 
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The boundary condition on spatial infinity is introduced 

n\ r =oo = n , (2.1) 

where r = ((a; 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 ) 1 ' 2 and no is a fixed vector, e.g. corresponding to 
the north pole 

n = (0,0,1). 

We shall consider mostly the time independent configurations, corresponding to a 
soliton at rest. The boundary condition (2.1) effectively compactifies the space R 3 , 
turning it into sphere § 3 , thus the stationary configurations realize the map 

n:S 3 ^S 2 , (2.2) 

which are known to be classified by Hopf invariant, sort of topological charge. 

In general, the density of topological charge is the zero component Jo of the 
current J M , which is conserved 

= 

independently of the equations of motion. Mathematically it is more natural to use 
the 3-form J dual to 1-form J* = J M d x^ and define the topological charge as an 
integral of J over space section 

Q= f J. 

In our case the 3-form J is constructed as follows. The pull-back of the volume 
2-form on S 2 via map (2.2) defines the closed 2-form on the space time 

H = ff^d/Ad/, 

where antisymmetric tensor H^ u is expressed via field configuration n(x) as follows 

Hfj, v = (<9 M n x d v n, n). (2.3) 
Here I use usual notations of vector analysis in 3-space. In fact H is exact 

H = dC 

and current 3-form is given by 

J = ]-H A C. 

An 

In more detail, we have the relations 

Hik = diCk — dkCi 

and 

1 / 3 

Q = ~ — / £ikj HikCjd x. 

For regular configurations Q gets integer values. This integer has a nice inter- 
pretation in the description of which I shall use the terminology of magnetostatic. 
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Tensor ifjfc can be interpreted as a field strength of the stationary magnetic 
field in Maxwell theory. The corresponding lines of force are defined via equations 

d_ _ 1 

ds Xl ~ 2 6ikj kjj 

where s is a local parameter along the line. It is easy to see that components of 
n(x) along these lines are constant 

Ts n(x) = 0, 

giving two "integrals of the motion" . In other words, the Maxwell lines of force are 
the preimages of points on § 2 under the map ( |2.2| ). Hopf invariant is the intersection 
number of any pair of such lines. 

All these facts are well known and can be found in textbooks (see e.g. Jl4|). 
However I decided to include them into my text to make it more selfcontained. 



3. The dynamical model 

I introduce the dynamical model by giving the relativistic action functional 



A = a J {d tl nyd 4 x + bj {H^Yd^x. 

In the usual convention of high-energy physics A is dimensionless, so the parameter 
a has dimension [length] ~ 2 and parameter b is dimensionless. Corresponding static 
energy E has the same form as A with space-time coordinates substituted by space 
variables only 

E = a I {d k n) 2 d 3 x + b [ (H ik ) 2 d 3 x. (3.4) 



and has proper dimension [length]" 1 . The structure of E is similar to that of 
Skyrme model, where the field variable having values in § 3 is used and corresponding 
topological charge is just a degree of map. 

Usual check based on the scale transformation is favorable for (3.4) in the same 
way as in Skyrme model. Indeed 

E = E2 + £4, 

where E2 and £4 are quadratic and quartic in derivatives of n correspondingly. 
Thus under scaling x — > Ax we have 

E2 — > XE2 , E4 — > — E4 
A 

and the virial theorem states that on the minimal configuration (if any) 

E2 = E4. 
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In terms of quantum theory Ei has a standard interpretation of the energy 
of nonlinear sigma-model whereas E4 is rather exotic. On the contrary in the 
magnetic interpretation, mentioned above, £4 is a natural term — it is just the 
Maxwell magnetic energy, whereas the nature of E2 is not that clear. However in 
what follows the presence of both E2 and E4 is crucial for the existence of solitons 
as the scaling argument already showed. 

This is confirmed also by a beautiful estimate, obtained in ||i5fl 

E>c\Q\ 3 /\ 

which shows that in the sectors with nonzero Q the minimum of energy is strongly 
positive. Thus the soliton solutions should be obtained by the minimizing of E with 
Q ^ fixed. 

Unfortunately until now there exists no proof of the compactness of the mini- 
mizing sequence in general case. For the case of axial symmetry uncour aging result 
are obtained in |l6| ] . So the main argument for the evidence of solitons in my model 
is based on the numerical work. 



4. Numerical work 

To find the numerical evidence of the existence of localized solitons it is not 
necessary to solve the nonlinear elliptic equation, obtained by the variational prin- 
ciple 

f£-* (4.5) 
on 

Instead one can introduce an auxiliary time s and consider the parabolic equation 



dn _ SE 
ds 8ft 

with initial value rimit 

n L=0 = j^init 



(4.6) 



bein g a configuration with the prescribed Hopf invariant. Of course to simulate 
( |4.6|) o n the computer one is to use some difference scheme. If for large s solution 



of (4J3) stabilizes it gives the solution of (4.5). In other words the soliton appears 
as an attractor for the evolution equation. 

There are of course many important practical details how to discretize equa- 
tion, how to take into account the normalization condition n 2 — 1 and how to 



choose the initial configuration n; n it- The main papers 10 and |ll|] use different 
prescription for all this, however quite satisfactorily the final results coincide. I refer 
to these papers for the details of calculations and proceed to describe the results. 

The iterative process was performed for the configuration with Q = 1, 2, ... 7. 
The results are as follows: for Q = 1 and Q = 2 the solutions are axial symmetric. 
The center line — the preimage of the point n — (0,0,-1) — is a circle. The 
surfaces = a, — 1 < a < 1 are toroidal and they are spanned by the lines of force 
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wrapping the torus once for Q — 1 and twice for Q = 2. In other words the soliton 
can be viewed as a filament of lines of force, closed and twisted once or twice. 

The solution for Q = 3 is similar but not axial symmetric any more, the 
corresponding "cable" is warped. For Q = 4 the soliton is a link of two twisted 
filaments. Especially beautiful case is Q = 7, the central line of the corresponding 
soliton is a trefoil knot. 

The file Jl7[ contains impressive moving pictures illustrating the convergence 
of the iterations. I plan to show these movies in my talk, but unfortunately can not 
do it in a written text. 

Thus the numerical work gives the compelling evidence of the existence of 
string-like solitons in my model. There remains an important mathematical chal- 
lenge to provide the rigorous existence theorem. Another interesting direction is 
to find some realistic applications of the model. Some progress in this direction is 
already obtained and I proceed to the description of it. 



5. The applications 

Nonlinear fields such as n(x) rarely enter the dynamical models directly. How- 
ever they can appear as a part of degrees of freedom in a suitable parameterization 
of the original fields. For example in condensed matter theory one uses the complex 
valued functions ip a (x), a = 1, .. . , N to describe the density amplitudes of Bose gas 
or the gap function of superconductor. The interaction supports the configurations, 
for which 

JV 
a=l 

is nonvanishing. In this case it is natural to use p as one of the independent variables 
and introduce new variables 

Xa = 1p a /p 

such that 

N 
a=l 

In this way the compact target (I use the slang of the string theory) § 2JV_1 
appears. 

When magnetic interaction is introduced the invariance with respect to the 
phase transformation 

fl> a (x) -^e lX(x ^ a (x) 
is invoked. This means that the target S 2JV_1 changes 

g2JV-l § 2JV-1 /C7(1) _ 

In particular for N = 2 we have 

S 3 /U(l) - § 2 
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and the field n(x) naturally appears. Quite satisfactorily the tensor Hn~ also emerges 
as a contribution to the magnetic field strength. 

Let us illustrate it in more detail. From the beginning we shall treat the 
stationary system, so no electric field will be used. 

The magnetic field is described in a usual way by means of the vector potential 
A k (x) and its interaction with -0-fields is introduced via covariant derivatives 

Vfe-0 = dki> + iA k ip. 

The energy density (of Landau-Ginsburg-Gross-Pitaevsky type) looks as follows 

2 i 

E = J2N^ a \ 2 + ^Fg + V(\^ a \), (5.9) 

a=l 

where 

F lk = d. t A k - d k Ai 

is the field strength of the magnetic field. The energy is invariant with respect to 
the gauge transformations 

ip a ^e iX ip a , A k ^A k -d k \. 

We shall make the change of the field variables so that only ga uge invariant ones 
will remain. For that observe that the first term in the RHS of (5J3) is a quadratic 
form in A 

2 2 

i v ^«i 2 = J2 i^i 2 + AkJk + p 2A *' 

a—l a—1 



where we use variable p from (|5.7| ) and introduce current 



2 



Jk = -i^i^ad^a - dk1pai>a)- 



a=l 

It is easy to check, that under the gauge transformations the current J k changes as 
follows 

J k ^J k + 2p 2 d k \ 
so that the sum ^ 

C k = A k + —-zJ k 

is gauge invariant. We shall use this variable instead of A k . Another gauge invariant 
combination is given by the quadratic form 

ft={xi,Xi)?(l 1 ), (5.10) 



\2 



where f = (ti, T2, T3) are Pauli matrices 



1\ fO -i\ I 



ri = U ' T2 = U ' T3 -\o -1 
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Normalization (5.S) for \a implies that ft is a real unit vector. In fact the map 



defined in (5.1C) is a standard Hopf map. Variable ft is manifestly gauge invariant 
and the set of variables (p, n, Ck) is our gauge invariant choice, substituting for the 
initial set (ip a ,Ak). The energy density can be explicitly expressed via p,n and Ck 
as follows 

E = {d k pf + p 2 ((d k n) 2 + C 2 ) + ^(d k d - d.Ck + H lk ) 2 + v(p, n 3 ). 



The most notable feature is the appearance of the tensor H^, defined in (2.3). 
The model, described in the main text, emerges if we put p — const and C = 0. 
Hopefully nontrivial p and C, at least confined to some range, do not spoil the 
soliton picture. This problem is under discussion now, see Jl|, p8[ . 

Let us stress, that the use of two fields ip a , a = 1,2 is most essential in this 
example. If N = 1 only variables p and C remain after the reduction, similar to 
just described. If N > 2 the CP(N — 1) field generalizing n has no topological 
characteristics. 

Another application, considered recently , deals with the parameterization 
for the SU(2) Yang-Mills field fi = 0, 1,2,3, a = 1,2,3. The Yang-Mills 

Lagrangian is invariant with respect to the nonabelian gauge transformations 

5A« = d^ a + f ahc A\e c . 

However in some treatments one reduces this invariance by the partial gauge fixing 
to the abelian one 

8B^ = ieB^ , SAl = d^e, 

where = A 1 ^ + iA 2 ^ is a complex vector field. I shall not discuss the reason for 
this reduction here and proceed assuming that it is done. Observe, that two vector 
fields A 1 ^, A 2 ^ in generic situation define a plane in Minkowski space and introduce 
an orthonormalized basis in this plane e", a = 1,2. 

Let e M = + ie 2 ^. The basis is defined up to rotation 



The fields B^ can be written in terms of this basis as 

and thus two complex valued fields and ip2 appear. The situation becomes quite 
similar to the previous example and indeed in [12] the complete parameterization of 
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the Yang-Mills variables is introduced with appearance of n-field and correspond- 
ing iJ-tensor. This is an indication that the Yang-Mills theory can have string-like 
excitations. However the situation is not that simple. The classical Yang-Mills the- 
ory is conformally invariant and has no dimensional parameters. Thus no hope for 
the localized regular classical solution exists. Nevertheless this complication could 
be lifted by quantum corrections. The famous "dimensional transmutation" , which 
leads to the appearance of dimensional parameter in quantum effective action, could 
favor the nonvanishing value of the corresponding p-variablc. All these considera- 
tions at the moment are rather speculative and need much more work to become 
reasonable. Personally I am quite impressed by this possibility and continue to work 
on it. 

6. Conclusions 

I think that the topic of my talk is quite instructive. It connects different do- 
mains in mathematics and mathematical physics: nonlinear PDE, elementary topol- 
ogy, quantum field theory, numerical methods. It illustrates the essential unity of 
mathematics, theoretical and applied. Finally it could lead to the realistic physical 
applications. For all these reasons I decided to present it to the ICM2002. 
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